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He II optical depth and UV escape fraction of galaxies 
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ABSTRACT 

We study the effect of H i ionizing photons escaping from the high-z galaxies on the 
Hen ionizing ultraviolet background (UVB) radiation. We show, while these photons 
do not directly interact with He n ions, they play an important role through radiative 
transport in modifying the shape of He n ionizing part of UVB spectrum. Within the 
observed range of UV escape from galaxies, we show the rapid increase in Hen Ly 
a effective optical depth at z ~ 2.7 can be naturally explained without resorting to 
pre-overlap era of He n reionization. Therefore, a well measured He n Ly a effective 
optical depth vs z relationship can be used to constrain the redshift evolution of UV 
escape from high-z galaxies. Our study also stresses the importance of including galaxy 
contribution even in the fluctuating UV background calculations. 
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1 INTRODUCTION 

The spectroscopy of z ^ 6 QSOs (Fan et al. 2006) and 
cosmic microwave background (CMB) polarization observa- 
tions (Larson et al. 2011) suggest that H I in intergalactic 
medium (IGM) got reionized at 6 ^ z ^ 12. The reionized 
IGM at z ^ 6 is believed to be in ionization equilibrium with 
UVB emanating from QSOs and galaxies (Haardt & Madau 
2012, HM12 hereafter). Measured temperature of the IGM 
at z ^ 3 (Becker et al. 2011), claimed excess temperature 
around z ~ 6 QSOs (Bolton et al. 2012), lack of substan- 
tial contribution to E > 54.4 eV photons from galaxies and 
redshift distribution of QSOs (Hopkins et al. 2007) all favor 
He II reionization occurring over 2.5 ^ z ^ 6.0. 

Direct measurements of He II Ly-a absorption from 
IGM is possible towards few UV bright high-z QSOs us- 
ing Hubble space telescope (for summary of observations 
see, Shull et al. 2010). The cosmic- variance limited available 
data suggest a rapid evolution of He II Ly-a effective optical 
depth (r ai Heii) and large fluctuation in the column density 
ratio of He II and H I (called as rf) over 2.7 ^ z ^ 3.0. The 
rapid evolution of the T a ,Heii is attributed to the comple- 
tion of He II reionization at this epoch (Furlanetto & Oh 
2008; McQuinn et al. 2009). The large fluctuation in r/ over 
small scales can be attributed to (i) small number of bright 
QSOs within a typical mean free path (A m f p ) of He II ioniz- 
ing photons (Fardal et al. 1998; Furlanetto 2009); (ii) large 
scatter in the QSO spectral index (Shull et al. 2004); (iii) 
presence of collisionally ionized gas (Muzahid et al. 2011) 
and (iv) small scale radiative transport effects (Maselli & 
Ferrara 2005). Most of the theoretical calculations of He II 
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optical depth consider only the QSO emissivity and radiative 
transport assuming the IGM gas in photoionization equilib- 
rium. While galaxies do not contribute directly to the He II 
ionizing radiation they can influence the ionization state of 
the IGM gas thereby affecting the He II optical depth. Here, 
we explore this using a cosmological radiative transfer code 
similar to HM12. We show that the A m fp is very sensitive to 
escape fraction (/esc) of H I ionizing photons from galaxies 
and the redshift evolution of r a; H C n can be used to con- 
strain / csc from high-z galaxies. Details of our calculations 
assuming (f2 m ,f2A,/0 = (0.3, 0.7, 0.7) cosmology is given in 
the following section. 



2 UVB CALCULATION 

We calculated the UV background spectrum contributed by 
QSOs and galaxies using the standard assumption that each 
volume element is an isotropic emitter and sink (see for ex- 
ample, Haardt & Madau 1996; Fardal et al. 1998; Faucher- 
Giguere et al. 2009). The angle and space averaged specific 
intensity J„ (in units of erg cm" 2 s" 1 Hz" 1 sr" 1 ) of diffuse 
UV background as seen by an observer at a redshift zo and 
frequency uq is given by (Haardt & Madau 1996): 
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where v = foil + z)/(l + Zo) is the frequency of emitted 
radiation at redshift z, e„(z) is the proper space averaged 
specific volume emissivity of radiating sources (QSOs and 
galaxies), 4^ is the proper line element in the Friedmann- 
Robertson- Walker cosmology and r e g is the effective optical 
depth which quantifies the attenuation of photons observed 
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at a frequency v§ while travelling through the IGM in be- 
tween z and zq. 

If we assume that the IGM clouds with neutral hydro- 
gen column density, JVhi, are Poisson-distributed along the 
line of sight, we can write T e a as (Paresce et al. 1980), 



Teff (^o, Zo, z) = I dz' 



d 2 N 
dNm dz' 



. (2) 



where, /(JVhi, z) = d 2 N / dNuidz, is the number of absorbers 
with JVhi per unit redshift and column density interval mea- 
sured at z. This is directly measured with QSO spectroscopy 
(see Petitjean et al. 1993). Assuming that absorbing clouds 
are made up of a pure H and He gas, the continuum optical 
depth through an individual cloud can be written as, 

t{v') = N m am(y') + NH elCHeltf') + AfHeIlCHeIl(^') • (3) 

Here, v' = 2Vq(1 + z')/(l + Zo) and N x and a x are the col- 
umn density and photoionization cross-section for a species 
x respectively. From QSO spectroscopy we know JVhcI is neg- 
ligible when A?hi ^ 10 17 ' 2 cm -2 . Even for Lyman limit sys- 
tems (Nm > 10 17 ' 2 cm -2 ), the ratio Nuei/Nm is small (see 
HM12) enough to neglect its contribution over the redshift 
range of interest in our study. Therefore, Eq. (3) becomes 



t{v) « Nm[<rai{v) + V (XHell(^)]- 



(4) 



In the absence of direct measurements of f(Nu c ii,z), a 
knowledge of 77 as a function of Nm together with f(Nm, z) 
allows us to calculate the contribution of He II to the con- 
tinuum optical depth. Under photoionization equilib- 
rium the JVhi and Nh c ii are related through the fol- 
lowing quadratic equation (Fardal et al. 1998; Faucher- 
Giguere et al. 2009; Haardt & Madau 2012), 

"He ^HI r 912,HI _ ^Hell T"228,HeII 

~A ft 1 A \ ~ T 228,HcII + —5 r {?) 

4n H (1 + A T 9 i2,Hl) (1 + B T 2 28,HcIl) 

where, t\, x is N x a x (X) for the species x. The values of A 
and B depend on the assumed relationship between the 
total hydrogen column density (N-jj) and electron den- 
sity (n e ). This relationship is obtained for a constant 
density slab of gas with a thickness of Jeans length un- 
der optically thin photoionization equilibrium (Schaye 
2001). Numerical simulations suggest that such a re- 
lationship is valid for log A^h^18 (see Rahmati et al. 
2012). We have considered a series of photoionization 
models of plane parallel slab having N^-n e relation- 
ship of Schaye (2001) illuminated by a powerlaw source 
using the photoionization code Cloudy (see Ferland 
et al. 1998). We confirm that the values A = 0.02 and 
B = 0.25 (as used by HM12) provide a good fit to the 
model predictions and adopt these values in our cal- 
culations. The quantity I x for x th species is defined as 
Ix = F x /n e a x (T), where a x (T) is the case-A recombi- 
nation coefficient and T x is photoionization rate of X th 
species as given by, 



dv 4ir a x {v) 



(6) 



Here v x is the ionization threshold frequency for the 
species x. We use T = 2 x 10 4 K and the form of n e and 
f(Nni, z) as given in HM12 in all our calculations. 
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Figure 1. UV background spectrum at z = 3 from our models for 
different values of k. The spectrum plotted in red is from HM12. 



2.1 Quasar and Galaxy emissivity 

The specific volume emissivity of radiating sources, 
62,(2:), is a sum of the emissivity from galaxies, e Ut a(z), 
and quasars, c v ^q{z). We have considered the paramet- 
ric form to the observed quasar co-moving emissivity at 
912A as given in HM12, 



E912,q(z) 

(1 + z) 3 
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24.6 
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4.68 



exp(-0.28z) 
exp{l.llz) + 26.3 



(7) 



in units of ergs Mpc" 3 s" 1 Hz" 1 together with the 
broken power law spectral energy distribution (SED), 
L v oc iv -044 for A >1300A and L u oc i/ -1 ' 57 for 
A < 1300 A (Vanden Berk et al. 2001; Telfer et al. 2002). 
To calculate the co-moving emissivity from galaxies we 
have taken the parametric form of star formation rate 
density (SFRD) used in HM12 : 



SFRD(z) 



6.9 x 10" 3 + 0.14(z/2.2) 
1 + (z/2.7) 4 - 1 



1.5 



M Q yr _1 Mpc~ 



( 8 ) 

The co-moving emissivity of galaxies (in units ergs 
Mpc" 3 s" 1 Hz" 1 ) is taken to be, 



(I + 2:) 3 



= C(z) x SFRD(z) x l u (z,Z) 



(9) 



Here, l v (z, Z) is the specific luminosity of a galaxy pro- 
duced for every solar mass of gas having metallicity Z 
being converted to stars. We get l v (z, Z) using the stel- 
lar population synthesis model 'Starburst99 V6.0.3' 
(Leitherer et al. 1999) for Z=0.001 and Salpeter initial 
mass function with stellar mass range 0.1 to 100 M Q . 
SED fitting studies, semi-analytic modelling of luminos- 
ity function and spatial clustering are consistent with 
star formation in a typical galaxy lasting for a few 100 
Myr (see discussions in Jose et al. 2012, and references 
there in). Therefore, for simplicity we assume that the 
star formation has lasted for more than 100 Myr. This 
assumption allows us to get a linear relationship be- 
tween star formation rate and luminosity. We do not 
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Figure 2. rj as a function of N//j for redshifts 2 (blue), 2.5 (black) 
and 3 (red). The dash, solid and dot-dash curves are for k = 0, 1 
and 2 respectively. 
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Figure 3. Hi Lya effective optical depth as a function of z for 
different values of b parameters. The observed points are from 
Becker et al. (2012). 



include recombination emissivity and resonant absorp- 
tion effects in our calculations. 

The factor C(z) is used to modify the SED to take 
care of the reddening and UV escape. For A >912A, 
C(z) is exp(—T u ) with t v being the frequency depen- 
dent dust optical depth calculated using extinction law 
of Calzetti et al. (2000) with R v = 3.1. The dust correc- 
tion may depend on the luminosity and z of individual 
galaxies (see for example Bouwens et al. 2012). How- 
ever, for simplicity we use a single value at all z. The 
dust optical depth is chosen to have a reddening correc- 
tion factor of 3 at 1500A. 

For wavelength range 228A< A <912A we take 
C(z) — f csc . To take into account the redshift evolution 
in the / C sc found by Inoue et al. (2006) we adopt the 
form of HM12 and use, 

C(z) = / esc = k [3.4 x 1(T 4 (1 + z) 3 - 4 ] . (10) 

In this case C{z) corresponds to the absolute escape 
fraction, / osc , defined as the ratio of escaping Lyman 
continuum (LyC) flux from a galaxy to the one which is 
intrinsically produced by the stars in it (Leitherer et al. 
1995) and k is a free parameter that allows us to change 
the values of / esc - For our fiducial star formation model, 
at k = 1 we get UV emissivity from galaxies similar to 
that of HM12. Models with fc = corresponds to spec- 
trum contributed by QSOs alone at A <912A. Note, 
no reddening correction is applied for A < 912A and 
we simply scale the unattenuated spectrum by / csc . In 
other words, we assume all the ionizing photons that are 
escaping through holes in the galaxy. We take C(z) = 
for A <228A because galaxies at the redshifts of our in- 
terest do not produce sufficient Hen ionizing photons 
as massive Population-Ill stars are rare. 



3 RESULTS AND DISCUSSION 

In Fig. 1, we present the UVB spectrum at z = 3 com- 
puted using our numerical code for three different k 



values. At z — 3, k = l and 2 corresponds to / csc of 
~4% and ~8% respectively. This is well within the in- 
ferred range for high-z galaxies and Ly-a emitters at 
z ~ 3 (Shapley et al. 2006; Iwata et al. 2009; Boutsia 
et al. 2011; Nestor et al. 2012). The range of r_i 2 , H i = 
0.4 — 1.3 predicted by our models is also consistent with 
the range allowed by the observations (see, Meiksin & 
White 2004; Bolton & Haehnelt 2007; Faucher-Giguere 
et al. 2008). For comparison we also plot the UVB spec- 
trum generated by HM12. It is clear that our model 
with k = 1 reproduces HM12 spectrum very well. The 
differences in Thi and THeii between the two codes is 
much less than the spread in these values due to allowed 
range in / csc . It is clear from Fig. 1 that increasing (de- 
creasing) k increases (decreases) Thi and decreases (in- 
creases) ThcII- As A $;228A emissivity is not affected 
by galaxies, the variation of Tneii with / esc can be at- 
tributed to the effect of / esc on He n opacity. Below we 
explore this in more detail. 



3.1 Escape fraction and r/ 

In Fig. 2 we have plotted r\ vs. Nhi at redshift 2, 2.5 and 
3 for different k values. Two trends are clearly evident 
for log A^hi ^ 17 : (i) for any given / CS c, V increases with 
increasing z and (ii) at any z, r\ increases with increas- 
ing / osc . The first trend can be understood in a simple 
way for the optically thin limit where r\ oc Fhi/FhcII- 
Both QSOs and galaxies contribute to Thi but to ThcII 
only QSOs alone contribute. Therefore, r] depends on 
how population of QSO and galaxy evolve with red- 
shift. For z ^ 2.5, population of QSOs declines rapidly 
(Ross et al. 2012) while that of galaxies remains almost 
same (Bouwens et al. 2011) helping n to increase. How- 
ever, at any z by increasing / GSC , we are increasing the 
galaxy contribution to Fhi which will further increase 
77. This explains the second trend we notice in Fig. 2. A 
simple physical picture we can draw is, a given Ngi is 
produced by integrating over larger column length when 
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Figure 4. Lya effective optical depth for Hell as a function of 
redshift for different f CS c and ft parameters (ft = 28 dash curve, b = 
30 solid curve, b = 32 dot dash curve) in km/s. Black diamonds 
are observations of Worseck et al. (2011). 



Tfji is high and therefore we get more -/Vneii- However, 
the implied ionization corrections for the range of / esc 
considered here still keeps fi;, within the allowed range 
from CMB observations. 



3.2 Ly a effective optical depth for Hi and Hell 

Now we estimate the Ly a effective optical depth of 
Hi (i.e., r ai m) and Hen (i.e., r ai H e ii) for a range of 
/ csc and compare with the available observations. The 
t Q..HI and To^HeH is given by (Paresce et al. 1980; Madau 
& Meiksin 1994), 



1 



dN m f(N m ,z)W n 



(11) 



where, X UiX is 1215.67A for Hi , 303.78A for Hen and 
W n is equivalent width of corresponding line expressed 
in wavelength units given by, 



-v<t>W\ 



(12) 



Here y = Nm for Hi , y = r/Nm for Hen and (j)(X) is 
Voigt profile function. In Fig. 3, we show our best fit to 
the observed T a ,m vs z relationship (Becker et al. 2012) 
obtained for b parameter values b = 30 ± 2 kms . 

In Fig. 4, we have plotted t q ,h c ii vs - z for different 
/osc values along with observations of Worseck et al. 
(2011). Here, we assume that the non-thermal motions 
dominate the He n line broadening and used the best fit 
values of b obtained for H I. In all cases r^Heii increases 
with increasing z. However the rate of increase depends 
on / esc . Interestingly, a sharp raising trend of r Q .Heii at 
z ~3, seen for k = l, almost provides a lower envelop 
to the observations. For k = 2, the model prediction al- 
most passes through most of the observed mean points 
at z < 2.7. This means that / e sc can not be much higher 
than 8% in this redshift range. To check whether the 
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Figure 5. Mean free path for He II ionizing photons as a function 
of z for different k values. Here, distance is in proper length scale. 



rapid rise in T a ^ c ii is a consequence of (1 + z) 3A evolu- 
tion assumed for / csc (Eq. 10), in Fig. 4 we also plot the 
Ta.Heii computed for redshift independent / csc (i.e 4% 
and 8%). Even in these cases, the z evolution of T a) Heii 
is more steeper than for k — 0. This observed rise and 
scatter in r aj H e ii is usually attributed to the pre-overlap 
era of He III bubbles assuming that the He II re- ionization 
completes around z — 2.7 (Dixon & Furlanetto 2009; Shull 
et al. 2010; Worseck et al. 2011). Our results suggest that the 
observed trend of T a! Heii with z can also be produced natu- 
rally by the radiative transfer effects associated with chang- 
ing /esc- Also a well measured lower envelop of r a! H e ii with 
z can be used to constrain the redshift evolution of / eS c. 

Next we calculate the A mp f predicted by our models. 
We define, 



A 



mfp 



dz 
dr cB 



(13) 



^ = / dNm/O/Vm, z) [1-e 



r(228A)1 



(14) 



where from Eq. (2) we write, 
dr c g 
dz Jo 

Here, r(228A) = N H i(cr HI 228A + ^ Hc ii,228a)- In Fi S- 5 the 
A m fp of Hell ionizing photons is plotted as a function of 
2 for different / osc . These curves are well approximated by 
Ampf(-z) = Ai x exp[— (z — 2)/Az] with best fit values of 
Ai = [150, 120, 90] Mpc and Az=[0.60, 0.41, 0.37] respec- 
tively for fc=[0,l,2]. This clearly demonstrates the rapid re- 
duction in A mp f with / ea c. This is a consequence of the fact 
that when we increase / csc , a given A^Heii is produced by a 
cloud with lower JVhi (i.e rj increases). As f(Nm,z) is pow- 
erlaw in iVm with a negative slope the A mtp reduces steeply. 
For our fiducial model with k = l, A mtp is 22 and 11 proper 
Mpc at z = 2.7 and 3.0 respectively. These are at least a 
factor 2 smaller than the corresponding values for fc = 0. 

Recently, Davies & Furlanetto (2012) have computed 
the mean evolution of r a! H e ii assuming diffuse emissivity 
of QSOs and allowing for fluctuating FhgII. This model also 
produces a rapidly evolving Tq^hcii with z without implicitly 
assuming He II reionization around z ~ 3 when a minimum 
A m fp of 35 comoving Mpc is used. This is already more than 
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A m fp we get at z — 3 for our fiducial model. Therefore, we 
can conclude that the galaxy contribution, that is usually 
ignored, will play an important role in fluctuation calcu- 
lations. We also speculate that the redshift at which the 
fluctuations beginning to dominate will also depend on / esc . 
Detailed investigations of high ionization species like N v 
and O VI absorption in QSO spectra at z ^ 2.5 and prox- 
imity effect analysis of He n Ly-a forest will provide more 
insights into the issues discussed here. We plan to address 
these carefully in near future. 



4 CONCLUSIONS 

We studied the effect of escape fraction of H I ionizing pho- 
tons from high- 2 galaxies on the UVB calculated by solving 
cosmological radiative transfer equation. We have demon- 
strated that the H I ionizing photons from galaxies play an 
important role in deciding the shape of He II ionizing part 
of UVB and affect the transmission of UVB through 
IGM. The main results are summarized below : (i) 
Hen ionizing part of UVB depends strongly on the 
value of /esc for z > 2.0. This is a consequence of the de- 
pendence of ratio n = Nn e u/Nui on / csc . (ii) Mean free 
path of Hen ionizing photons decreases rapidly with 
increasing / csc . Suggesting the fluctuations in 77 can be 
much higher compared to the one computed ignoring 
the galaxy contributions, (iii) We show for the range 
of / esc allowed by the observations the rapid increase 
in Hen Lya effective optical depth, recently observed 
at z ~ 2.7, can be naturally explained without invoking 
pre-overlapping era of He 11 re- ionization around this 
redshift. Thus one will be able to place strong con- 
straints on /esc as a function of z using effective optical 
depth evolution of Hen . We also show the observa- 
tions at z ^ 2.7 are consistent with / csc being not much 
higher than 8%. 
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